Development of software for human muscle force estimation

Muscle force estimation (MFE) has become more and more important in exploring principles of pathological movement, studying functions of artificial muscles, making surgery plan for artificial joint replacement, improving the biomechanical effects of treatments and so on. At present, existing software are complex for professionals, so we have developed a new software named as concise MFE (CMFE). CMFE which provides us a platform to analyse muscle force in various actions includes two MFE methods (static optimisation method and electromyographic-based method). Common features between these two methods have been found and used to improve CMFE. A case studying the major muscles of lower limb of a healthy subject walking at normal speed has been presented. The results are well explained from the effect of the motion produced by muscles during movement. The development of this software can improve the accuracy of the motion simulations and can provide a more extensive and deeper insight in to muscle study.     

Surface EMG-force modelling for the biceps brachii and its experimental evaluation during isometric isotonic contractions

The aim of this study was to evaluate a surface electromyography (sEMG) signal and force model for the biceps brachii muscle during isotonic isometric contractions for an experimental set-up as well as for a simulation. The proposed model includes a new rate coding scheme and a new analytical formulation of the muscle force generation. The proposed rate coding scheme supposes varying minimum and peak firing frequencies according to motor unit (MU) type (I or II). Practically, the proposed analytical mechanogram allows us to tune the force contribution of each active MU according to its type and instantaneous firing rate. A subsequent sensitivity analysis using a Monte Carlo simulation allows deducing optimised input parameter ranges that guarantee a realistic behaviour of the proposed model according to two existing criteria and an additional one. In fact, this proposed new criterion evaluates the force generation efficiency according to neural intent. Experiments and simulations, at varying force levels and using the optimised parameter ranges, were performed to evaluate the proposed model. As a result, our study showed that the proposed sEMG–force modelling can emulate the biceps brachii behaviour during isotonic isometric contractions.     

Analysis of tendinous actuation in balancing the maximal fingertip force for normal and abnormal forefinger system

This work displayed the force capabilities of the musculoskeletal system of the forefinger under external loading. Different states of normal and pathological fingers are studied. We evaluated the impact of losing musculo-tendon unit strength capacities in terms of maximal output fingertip force and tendon tensions distribution. A biomechanical model for a static force analysis is developed through anatomical and kinematic studies. An optimisation approach is then used to determine tendon tension distribution when performing an isometric task. Furthermore, pathological fingers with common cases of injured flexors and extensors are analysed. The method of simulation for forefinger abnormities is described. Furthermore, the simulation results are interpreted.     

Molecular dynamics simulations of a branched tetradecasaccharide substrate in the active site of a xyloglucan endo-transglycosylase

Molecular dynamics simulations of the tetradecasaccharide XXXGXXXG in complex with the hybrid aspen xyloglucan endo-transglycosylase PttXET16-34 have been performed and analysed with respect to structure, dynamics, flexibility and ligand interactions. Notably, the charge state of the so-called ‘helper residue’ aspartate 87 (Asp87), which lies between the catalytic nucleophile [glutamate 85 (Glu85)] and general acid/base (Glu89) residues on the same beta strand, had a significant effect on PttXET16-34 active site structure. When Asp87 was deprotonated, electrostatic repulsion forced the nucleophile away from C1 of the sugar ring in subsite ? 1 and the proton–donating ability of Glu89 was also weakened due to the formation of a hydrogen bond with Asp87, whereas the protonation of Asp87 resulted in the formation of a hydrogen bond with the catalytic nucleophile and correct positioning of the catalytic machinery. The results suggest that catalysis in glycoside hydrolase family 16, and by extension clan GH-B enzymes, is optimal when the catalytic nucleophile is deprotonated for nucleophilic attack on the substrate, whereas the ‘helper residue’ and general acid/base residue are both in their conjugate acid forms to align the nucleophile and deliver a proton to the departing sugar, respectively.     

Controlling the secondary-structure-forming tendencies of proteins by a backbone torsion-energy term

We examined a new backbone torsion-energy term proposed by us in the force field for protein systems. This torsion-energy term is represented by a double Fourier series in two variables, namely the backbone dihedral angles φ and ψ. It gives a natural representation of the torsion energy in the Ramachandran space in the sense that any two-dimensional energy surface periodic in both φ and ψ can be expanded by the double Fourier series. We can then easily control secondary-structure-forming tendencies by modifying the torsion-energy surface. For instance, we can increase or decrease the α-helix-forming-tendencies by lowering or raising the torsion-energy surface in the α-helix region and likewise increase or decrease the β-sheet-forming tendencies by lowering or raising the surface in the β-sheet region in the Ramachandran space. We applied this torsion-energy modification method to six force fields, AMBER parm94, AMBER parm96, AMBER parm99, CHARMM27, OPLS-AA and OPLS-AA/L, and demonstrated that our modifications of the torsion-energy terms resulted in the expected changes of secondary-structure-forming tendencies by performing folding simulations of α-helical and β-hairpin peptides.     

Nanoscale mapping and organization analysis of target proteins on cancer cells from B-cell lymphoma patients

CD20, a membrane protein highly expressed on most B-cell lymphomas, is an effective target demonstrated in clinical practice for treating B-cell non-Hodgkin's lymphoma (NHL). Rituximab is a monoclonal antibody against CD20. In this work, we applied atomic force microscopy (AFM) to map the nanoscale distribution of CD20 molecules on the surface of cancer cells from clinical B-cell NHL patients under the assistance of ROR1 fluorescence recognition (ROR1 is a specific cell surface marker exclusively expressed on cancer cells). First, the ROR1 fluorescence labeling experiments showed that ROR1 was expressed on cancer cells from B-cell lymphoma patients, but not on normal cells from healthy volunteers. Next, under the guidance of ROR1 fluorescence, the rituximab-conjugated AFM tips were moved to cancer cells to image the cellular morphologies and detect the CD20-rituximab interactions on the cell surfaces. The distribution maps of CD20 on cancer cells were constructed by obtaining arrays of (16×16) force curves in local areas (500×500 nm²) on the cell surfaces. The experimental results provide a new approach to directly investigate the nanoscale distribution of target protein on single clinical cancer cells.     

Direct force measurement of single DNA–peptide interactions using atomic force microscopy

The selective interactions between DNA and miniature (39 residues) engineered peptide were directly measured at the single‐molecule level by using atomic force microscopy. This peptide (p007) contains an α‐helical recognition site similar to leucine zipper GCN4 and specifically recognizes the ATGAC sequence in the DNA with nanomolar affinity. The average rupture force was 42.1 pN, which is similar to the unbinding forces of the digoxigenin–antidigoxigenin complex, one of the strongest interactions in biological systems. The single linear fit of the rupture forces versus the logarithm of pulling rates showed a single energy barrier with a transition state located at 0.74 nm from the bound state. The smaller k_off compared with that of other similar systems was presumably due to the increased stability of the helical structure by putative folding residues in p007. This strong sequence‐specific DNA–peptide interaction has a potential to be utilized to prepare well‐defined mechanically stable DNA–protein hybrid nanostructures. Copyright © 2013 John Wiley & Sons, Ltd.